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Abstract 



We compute isospin-violating meson-nucleon coupling constants and their 
consequent charge-symmetry-breaking nucleon-nucleon potentials. The cou- 
plings result from evaluating matrix elements of quark currents between nu- 
\ cleon states in a nonrelativistic constituent quark model; the isospin violations 

arise from the difference in the up and down constituent quark masses. We 
find, in particular, that isospin violation in the omega-meson-nucleon vertex 
dominates the class IV CSB potential obtained from these considerations. We 
evaluate the resulting spin-singlet-triplet mixing angles, the quantities ger- 
mane to the difference of neutron and proton analyzing powers measured in 
elastic n — p scattering, and find them commensurate to those computed orig- 
inally using the on-shell value of the p-oj mixing amplitude. The use of the 
on-shell p-oo mixing amplitude at = has been called into question; rather, 
the amplitude is zero in a wide class of models. Our model possesses no con- 
tribution from p-uj mixing at (7^ = 0, and we find that omega- meson exchange 
suffices to explain the measured n — p analyzing power difference at 183 MeV. 
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I. INTRODUCTION 



The suggestion of Goldman, Henderson, and Thomas [jl| that the contribution of p-uj 
mixing to charge-symmetry-breaking (CSB) observables is suppressed in the low momentum 
transfer regime has opened the search for new sources of isospin violation. Since then, 
many calculations, using a variety of models, have confirmed the suppression of the p- 
uj mixing amplitude at small spacelike momenta ^-0. Indeed, it has been shown that 
the p-uj mixing amplitude is zero at = in all models with vector mesons coupled 
to conserved currents 0. Yet, in Refs. 0-0] no alternate mechanisms to p-uj mixing are 
proposed. The phenomenological impact of this gap must be emphasized: the CSB potential 
from p-uj mixing — with the mixing amplitude fixed at the omega-meson point — accounts 
for some 40% of the difference between the neutron and proton analyzing powers (AA) 
measured in elastic n — p scattering at 183 MeV [||. Without this contribution the previous 
agreement between theory and experiment would be upset . Although the suppression 



of the mixing amplitude continues to be controversial [jl3|,|n|, sources of additional isospin 



violation are interesting in their own right and deserve examination. Indeed, the aim of 
the present paper is to show that a recently proposed CSB mechanism — based on isospin- 
violating meson-nucleon coupling constants [|l5l — is sufficient to restore the agreement 
with experiment. Specifically, we examine the effect of these new sources of CSB on the 
spin-singlet-triplet mixing angles; these are the fundamental dynamical quantities driving 

Most theoretical efforts devoted to understanding CSB observables use a nucleon-nucleon 
(NN) interaction constrained by two-nucleon data P, p!0| , |T6i . In such a picture, isospin vio- 
lations arise from electromagnetic effects and hadronic mass differences. Sources of CSB can 
be classified in terms of three distinct contributions: (i) isovector-isoscalar mixing in the me- 
son propagator, (ii) isospin-breaking in the nucleon wave function, and (iii) isospin breaking 
in the meson-nucleon and photon-nucleon vertices. Rho-omega mixing, the proton-neutron 
mass difference, and the difference between the electric charge of the proton and the neutron 
are typical examples of (i), (ii), and (iii), respectively. The existence of isovector-isoscalar 
mixing, such as tt-t] and p-u mixing, is well established. For example, p-u mixing has 
been observed experimentally in e~^e~ tt+tt" measurements at the cj— meson production 
point [|T7|. However, the suggested suppression of the mixing amplitudes at small spacelike 
momenta lessens their impact on CSB observables. This is, in part, why we consider other 
sources of isospin violation in this paper. 

Isospin breaking in the nucleon wave function in a hadronic model is driven by the 
neutron-proton mass difference. Indeed, it is through this mechanism that charged-pion 
exchange dominates [0,0 the class IV potential |^ at moderate momentum transfers. 
Isospin breaking in the nucleon wave function can also arise in a quark model picture from 
the mixing of the nucleon to \ J'^ = 1/2"'"; T = 3/2) baryon states |]2T|. While undoubtedly 
nonzero, one expects the T = 3/2 components of the nucleon to be small due to the large 
mass difference between the nucleon and the A(1910) — the first P31 baryon. In contrast, 
the p-u mass difference is a mere 12 MeV. Thus, we turn to the meson-nucleon coupling 
constants as the possible sources of isospin violation demanded by data. 

While there have been calculations of isospin-violating meson-nucleon coupling con- 



stants [11,22-25], their impact on class IV CSB observables has only recently been con- 
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sidered Here, as earlier [0, we adopt a nonrelativistic quark model to calculate 

isospin breaking in the meson-nucleon coupling constants. In the model the coupling con- 
stants emerge from evaluating matrix elements of quark currents of the appropriate Lorentz 
and flavor structure between nucleon states. The isospin violations arise from the up-down 
quark mass difference. Radiative corrections to the vertices have also been evaluated and 
are found to be small [p6l , p7| . Here we study the phenomenological impact of Ref. JTSf in 
greater detail. In order to do this, we estimate the dependence of the isospin-breaking 
found in the vertices at = 0. 

We have organized the paper as follows. In Sec. |T| the model is introduced, and isospin- 
violating meson-nucleon coupling constants are computed. We show that in the = limit 
the couplings depend merely on the spin-flavor structure of the nucleon wave function; they 
are insensitive to the spatial components of the nucleon wave function. In Sec. |IT1| we use 
these findings to compute the resulting CSB potentials. In particular, we obtain a large 
contribution from omega-meson exchange to the class IV potential. We quantify the impact 
of isospin- violation in the NNuo vertex by computing the resulting spin-singlet-triplet mixing 
angles — these are the basic building blocks of AA. These results are presented in Sec. 0. 
Finally, we discuss the impact of our work in Sec. 0. 



II. ISOSPIN- VIOLATING MESON-NUCLEON COUPLING CONSTANTS 

We are interested in computing the coupling of an on-shell nucleon to the neutral mesons 
uj, p°, 77°, and (J. The off-shell vertices could engender additional isospin breaking, but our 
primary focus is on the NN system, so that we will not consider these effects further. 
The exchanged mesons couple to nucleon currents of the appropriate Lorentz character, 
and the meson-nucleon coupling constants emerge from evaluating the matrix elements of 
these currents in the quark model. The difference in the up and down constituent quark 
masses thus gives rise to isospin-violating meson-nucleon coupling constants. At = these 
couplings are determined from the spin and flavor structure of the nucleon wave function 
alone. In contrast, the couplings at = of the nucleon to the charged mesons are sensitive 
to the quark momentum distribution as well, and are, therefore, more model dependent []1T 



We shall consider the neutral- vector- meson-nucleon vertices first, as they are relevant to the 
A A measurement. The most general form for these on-shell iViV- meson vertex functions, 
consistent with Lorentz covariance and parity invariance, are 

- ^wAj;^. = [dNr + ^/^^"^^4m^) ' ^^^^ 
-'s^nA'^m, = -'3^N, [9%!' + ^/^^'^^4m^) ' ^^^^ 

where Qj^j^^^ (a = i^,p) are the isospin-averaged, phenomenological, meson-nucleon cou- 
pling constants, determined from fits to the NN phase shifts and to the properties of the 
deuteron p8|j29[] , and M^r is the nucleon mass. We compute the couplings, that is, g%{q^) 
and f^{q'^), by assuming that the NNa vertex functions can be related to the matrix el- 
ements of quark currents of the appropriate Lorentz and flavor structure between nucleon 
states, computed in the nonrelativistic quark model. Thus, 
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{N{p', s')\J^-'^\N{p, s)) = U{p', s')A%^U{p, s) , (2a) 
{N{p', s')\J^-'-\N{p, s)) = U{p\ s')K^^^U{v. s) . (2b) 

Here U{p, s) denotes a on-shell nucleoli spinor of mass Mat, momentum p, and spin s. We 
shall focus on the couplings at = 0, where q = p' — p, as the nonrelativistic quark model 
is best suited to an estimate in the static limit. The quark currents J'^'^ are 

J^'+ = \uYu + \d^''d, (3a) 
jm;- = ^^f^u - dYd . (3b) 

It is the quark vector current which is appropriate to the vector-meson-nucleon vertex; the 
second superscript (±) denotes its symmetry under the u ^ d flavor transformation. Note 
that the constituent quarks are assumed to be elementary: no quark form factors have been 
introduced. The isoscalar vector quark charge is 1/3, whereas the isovector vector quark 
charge is +1 for the up quark and —1 for the down quark. The charge assignments are made 
such that g"^ = 1, = 1, and g^ = —1, at = 0. 

Our model stems from the notion of vector dominance |^0|- Vector dominance presumes 
that a photon's interaction with a nucleon is mediated by the rho — or omega — meson. 
Here we argue that the coupling of the vector mesons themselves to the nucleon can be 
determined via matrix elements of the appropriate isospin components of the quark vector 
current. Our model does not predict the isosp in-conserving coupling constants g^^j^^', these 
must be extracted from phenomenological fits to two-nucleon data. However, the isospin- 
violating pieces, as well as the tensor-to-vector ratio, can be calculated within the model. 
Note that the vector dominance nature of our model implies that the quarks couple to 
conserved currents. We estimate the resulting coupling constants using the nonrelativistic 
quark model (NRQM); this is an additional assumption. 

The couplings g'^ and are functions of the meson four-momentum g^, though we shall 
focus on the couplings at = 0. In this limit the couplings are insensitive to the spatial 
component of the nucleon wave function; they follow directly from its spin and flavor content 
alone. In the SU{6) limit ||31|| , 

IpT) = ^(2kT dl) ~\uUl d^) - \ul d^) + 



2\u] di nt) - |ni dt '"!) - l'«T d^ui) + (4) 
2\diu^u^) - \d^ u] ui) -\d] uiu])' 



The neutron spin-up wave function, \n^), is obtained by exchanging the up and down quarks 
in the expression for \p^). The isospin violations arise from the difference in the up and 
down constituent quark masses. The couplings constants are obtained from computing the 
matrix elements found in the nonrelativistic reduction of Eq. (^ in the quark model; i.e., 



gN = T.9t{N]\l\N]) ■ ^^ = T.^^t{NU^m^), (5a) 



3 /'^p I 3 



g'N = T.9r{N]\i\N^) ■ %^ = E/^r(iVT|c^-|iVT)- (5b) 
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Note that we have introduced the quark magnetic moment fxf = gf/2mi, with the charges 
given in Table |. In the following presentation we discuss only the coupling of the nucleon 

to the meson, as an illustrative example. Our results, collected in Table 0, include the 

couplings to the other mesons as well. 

The vector coupling of the u;— meson to the nucleon is determined by simply counting 

the quark charges: 



+ 



(6a) 
(6b) 



The tensor coupling, in contrast, depends on the spin structure of the nucleon wave func- 
tion: 
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It is useful to construct isoscalar and isovector combinations at the nucleon level; i.e., 



where 



9n 

9o+9irz-- 



9n^{l-rz)^g^ + grrz 



r,) = fi^ + /x^r. 



1, 
1 
6m 



5 Am 

1 H n 

6 m 



i9o+fo) , i9i+fr) 



2M 



2M 



T, 



(7a) 
(7b) 

(8) 
(9) 

(10) 
(11) 



with 



M = -(M„ + Mp) 



m 



-(mrf + m„) 



Am = {nid 



m„ 



(12) 



The expression in Eq. ([TT|) is given to leading order in Am/m only. Note that isospin 
breaking in the / and g couplings is realized in the / alone and that the breaking in the uj 
tensor coupling is isovector in character. The uo — and p — vector couplings are isospin- 
conserving. The isospin-breaking in our model is connected to that of the electromagnetic 
form factors through our assumption of vector dominance; charge conservation protects 
the charge form factor from isospin-breaking at zero momentum transfer The tensor 

coupling is explicitly sensitive to the quark mass, as seen in Eq. (^, and the isospin-breaking 
corrections are generated by the up-down mass difference. In the constituent quark model 



Am > |3^; the up quark, which is lighter, has a larger magnetic moment than the down 
quark. Henceforth we shall adopt the choice M/3m = 1 in reporting the coupling constants. 
Our results are summarized in Table |1[ 

We now consider the isospin-conserving results. We find for the tensor-to-vector ratio 
that 
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These results are qualitatively consistent with the /n/qn ratios which emerge from phe- 
nomenological fits to the A^A^ interaction p8| , |29[] — recall that the Bonn B potential pa- 
rameters |2^, for example, are f%/g% = ^ and /n/Qn — This consistency is intimately 



connected to the NRQM's ability to describe the nucleon magnetic moments and to our 
assumption of vector dominance. In the NRQM, with M = 3m, the anomalous magnetic 
moment is purely isovector: Kp = 2 and Kn = —2. Note that = 1.79 and k'^^ = —1.91. 
These successes gives us confidence in using our model to compute the isospin-violating 
corrections to these coupling constants. 

For completeness, we shall now consider isospin breaking in the NNa and NNn^ vertices 
as well. We exclude the NNi] vertex from this discussion because the g^^^^ coupling constant 



is poorly constrained by A^A^ data The appropriate vertex functions are 

^9NNa^%N. = '9nn. (9%) 1 , (14a) 

NNtt 



9nnAnn.= 9nnJ9n)i'- (14b) 



We have assumed pseudoscalar, rather than pseudovector, coupling for the pion in order to 
be consistent with earlier calculations of charge-symmetry breaking P, [T0i[T6i|T8i|T9| . In our 



model, we connect the vertex functions to matrix elements of quark currents, so that 

{Nip', s')\J^\Nip, s)) = Uip', s')A^^^Uip, s) , (15a) 
{N{p', s')\j'\N{p, s)) = U{p', s')Kl^U{p, s) , (15b) 

where 

J^{q) = - uu -\ — dd , (16a) 
3 3 

J^(q) = -u-f^u--dj^d . (16b) 
5 5 

The charges = 1/3, = 1/5, and g^ = —1/5 have been chosen such that g'^ = 1, g^ = 1, 
and g^ = —1 when Am = 0. Evaluating the nonrelativistic reduction of Eqs. ([15|a) and 
(P^) in the quark model, we find 

g% = j2gm^ m^) , (17a) 



9n 



i=l 
3 



Y^^,UNUcTnN]) , (17b) 



i=i 

where we have defined /if = gf/2mi. From Eq. (p!7|a) we see that the sigma meson generates 
merely a spin-independent coupling to the nucleon in the nonrelativistic limit, so that there 
is no isospin breaking in the NNa vertex and no contribution from sigma exchange to the 
CSB potential. Thus, we will not consider sigma exchange further. However, the quark 
mass dependence contained in in Eq. (p]7|b) implies that the isospin-breaking in the pion 
case is finite. The breaking to 0{Am/m) is indicated in Table ||. Note, however, that the 
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computed breaking at = depends on the nature of the assumed pion-nucleon couphng. 
If we had chosen pseudovector couphng, rather, then no isospin breaking would resuh. The 
pseudovector current contains no quark mass dependent pieces in the nonrelativistic hmit. 
Thus, our prediction in the 7r° case is decidedly more model dependent than in the p° and uj 
channels. Moreover, in the latter case, the compatibility of the computed tensor-to-vector 
coupling constant ratios with the Bonn potential indicates that vector dominance, which 
we assume, has some phenomenological support. Note that in the case, there is no 
such independent support of our "pseudoscalar dominance" assumption. This concludes our 
discussion of isospin breaking in the A^A^-meson vertices. 



III. CHARGE SYMMETRY BREAKING POTENTIALS 



We shall now compute the CSB potentials which arise from the isospin- violating couplings 
computed in the previous section and tabulated in Table 0. In an one-boson exchange 
approximation, presuming the form of the isospin breaking found in the = results, we 
obtain the following CSB potentials for u, and 7r° exchange, respectively: 



v; 



CSB 



CSB 



r^(l)7,(2)r,(l)+7^(l)r,(2)r,(2) 



^c'sB = nsBr^(i)r,(2) 



^.(2) 



+ ^SB 



r^(l)7;.(2)r,(2)+7^(l)r^(2)r,(i: 



VcsB = V^Jsb7'(1)7'(2) r,(l) + r,(2) 



where = ia^^{p' — p)yj1M and we have defined 

,2 \ 
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NNoj 



fi9o 



^NNp 



/o7f 



g2 _ ^2 



yoi/i 



q2 _ ^2 



(18a) 
(18b) 
(18c) 



(19a) 
(19b) 
(19c) 
(19d) 



Isospin breaking in the meson-nucleon vertices give rise to the above CSB potentials, as 
per Eqs. ( |19a| )-( [TI?n[ ). The isosp in-conserving tensor coupling is nonzero in the case of the p 
vertex, so that an additional potential of strength ^'csb(?) arises. These contributions have 
been considered only recently ||15| . Yet the potentials of Eqs. ( |18a|) and ( |18c|) are identical 
in form to those generated by p-u and tt-t] mixing, respectively. That is. 



^CSB 



^CSB 



r'^(l)7^(2)r,(l)+7^(l)r^(2)r,(2) 



VcSb = ^csW(l)7'(2) r.(l) + r.(2) 



(20a) 
(20b) 



where 
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^c^sb(?) = , ,A p\H\uj) , (21a) 

Note that in Eq. (pia|) we introduce /^^v^, the phenomenological tensor couphng of the 
Bonn model [^. Rather than performing a nonrelativistic reduction of the potentials in 
Eqs. ( |18a|) -( p!8cD and Eqs. ( P0a|) - (p0E| ), we simply classify the former as, either, "po;— like" or 
"vrry— like" potentials. The effect of these new isospin- violating potentials on CSB observables 
can then be readily elucidated. For example, the contribution from omega-meson exchange 
is identical in structure to that from p-uj mixing and thus contributes as well to AA in elastic 
n — p scattering. Indeed, we now show that the contribution from omega-meson exchange 
is comparable in magnitude and identical in sign to the one obtained from p-uj mixing — if 
the mixing amplitude is fixed at its on-shell value. 



1. One-boson exchange potentials of the p-uj kind 



Potentials of the form given in Eq. ( |20aD give rise to class III and class IV CSB potentials. 



They are generated by the interference between the isospin-conserving vector coupling and 



the isospin-violating tensor coupling; note, for example, Eq. (|18a|) and the second term in 
Eq. ( |18b|) . Unlike the case of the omega, the isospin structure of rho exchange is not identical 
to that of p-uj mixing; they are related by exchanging ^"^(l) ^ Tz{'^)- Thus, rho exchange 
contributes to the class IV p-uo mixing potential with a sign opposite to that of the omega. 
No sign changes are necessary when computing its vrr]— like or class III puo—Wke contribution. 
Note that the contribution from rho exchange is small relative to that from the omega — 
this emerges despite the larger isospin-violating coupling associated with the rho vertex (see 
Table |Tp. The vector NN p coupling is simply small relative to that of the omega; in the 
Bonn potential 9%^^/ Q^^p ~ 27 |^. The relative importance of the various contributions 
can be estimated by computing the CSB potentials at = 0. Recall that in this limit 
the isospin-violating couplings are insensitive to the quark momentum distribution; they 
depend only on the spin-flavor symmetry of the wave function. Using the Bonn B potential 
parameters of Table |T| and a value for the quark- mass difference of Am = 4.1 MeV ||32|| , 



we obtain the following results at = 0: 

(22a) 
(22b) 
(22c) 

Several remarks are in order. First, the p-u mixing amplitude, if modeled via fermion 
loops necessarily vanishes at = in our model. Our model assumes vector domi- 
nance, so that the vector-meson-nucleon vertices are determined by the appropriate isospin 





= 0) = 


9' 

mi 


i 2.49 GeV-^ , 




= 0) = 


^NNp rp p ^ 

2 /oyi ~ 
m^ 


0.18 GeV-^ . 




= 0) = 


f NNp^NNu 

m^m^ 
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components of the quark electromagnetic current. Thus, the vector mesons couple to cur- 
rents that are conserved at the nucleon level, so that the above result follows @J^. At the 
= point, the charge-symmetry violation in our model comes purely from the vertex 
contributions. Note that the rho meson contribution to the latter is, indeed, small. It rep- 
resents merely a 7% correction to the contribution from one-omega exchange. Second, the 
strength of the CSB potentials generated from omega exchange is comparable in magnitude 
to those obtained from p-uj mixing if the on-shell value of the mixing amplitude is assumed, 
{p\H\uj)\q2=rni = —4520 ± 600 MeV^ ||3^. Note, moreover, that the uj and on-shell p-uj 



mixing contributions are identical in sign. Specifically, 

VtsB^e = 0) = _ W^A^^ ^^|ff|^^ ^ 2.07 GeV-2 , (23) 

A CSB potential of this magnitude is needed for a successful description of AA at 183 
MeV 1^. Summing our omega exchange contribution to the CSB potential and that from 
on-shell p-u mixing is not only internally inconsistent but also gives a final potential which 
is too large to fit the data — see Sec. Our results suggest that a class IV potential 
of the appropriate size is generated by isospin-violations in the NNu vertex, together with 
small corrections from rho exchange and off-shell p-u mixing. This is the central result of 
our paper. 



2. One-boson exchange potentials of the n-ri kind 



Potentials of the form given in Eq. (pOb|) generate class III CSB potentials exclusively. 
The Lorentz structure of the first term in Eq. ( |18b| ) differs from that of the tt-t] mixing and 
one-pion exchange potentials, so that it is convenient to perform a nonrelativistic reduction 
of all three contributions, i.e.. 



^csB = ■~^csb(5') ^j^M^) 
%SB = -^csb{q) ^j^M^) 

K^SB = -^CSb(9) (yi^J^^^ 



CTi ■cr2 + Si2(q)J [r,(l) + r,(2) 
CTi ■0-2 + 5i2(q)] [r,(l) + r,(2) 
2o-i-cr2-5i2(q)l k(l) + r,(2) 



(24a) 
(24b) 
(24c) 



where we have introduced the tensor operator 5*12 (q) = [3(cri ■ q)(cr2 ■ q) — cti ■ (T2]. We 
estimate the relative size of these contributions by evaluating them at = 0, noting table |T|: 



v^JsB(r = 0) 

V^'sb{q' = 0) 



2 - - ^^gj^g'l ^ 36.83 GeV-^ , 



-fS ^ 1.07 GeV" 



(25a) 
(25b) 



The one-pion exchange contribution dominates that of the rho; this is driven by the large 
p-n mass difference — recall m?p/m\ ^ 30. Note that the inclusion of the rho meson leads 
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to a reduction of the tensor and an enhancement of the spin-spin components of the pion- 
exchange potentiaL Unhke the vector meson case, we cannot readily compute the n-i] mixing 
amphtude at = in our model. That is, in the pion case, there is no conserved current, 
so that the = mixing can be nonzero p5H37||. Nevertheless, we can compare the results 
of Eq. (p5D with the "usual" ir-rj mixing potential: 

Vc^q' = 0) = -^-^^^^^{n\H\v) , , ^ 52.01 GeV"^ , (26) 

where we have input the TT-rj mixing matrix element evaluated at its on-shell point, 
((K\H\rj)\q2^^2 = —4200 MeV^ ||3^. The contribution from one-pion exchange is compa- 



rable to that from n-rj mixing. The n-rj mixing potential may seem slightly larger, but the 
NNrj coupling is ill-determined from two-nucleon data. Indeed, it is believed that the Bonn 
potential overestimates it — a current analysis based on ry-photoproduction data suggests 
couplings as low as gj^j^^/ATi < 0.5 (see also Ref P5[]). 



The CSB potentials from one-pion exchange have been computed previously in a nucleon 
model |T3]. Here the neutron-proton mass difference, AM, generates the breaking. In the 



specific case of the class III contribution coming from neutral pion exchange, the scale of the 
breaking is set by AM/2M. Thus, the isospin breaking in the quark model is substantially 
larger than in the nucleon model, i.e., 

3 Am\ / /I AM\ , , 

so that any CSB observable receiving an important contribution from n-r] mixing will also be 
affected by the exchange of neutral pions. The breaking we calculate in the NNn^ vertex is 
identical to the result of Mitra and Ross |TT]J2^ . Note that the exchange of charged pions — 
and rhos — generates a class IV potential which is important in the analysis of AA |]10|,|ll|. 
In the charged meson case, however, the relation between the isospin-violating couplings 
in the two models is not simple: it depends on the quark momentum distribution. Yet, 
under reasonable assumptions, both sorts of models seem to generate class IV potentials of 



comparable strength [11 



IV. RESULTS 

In this section we compute the CSB potentials for a range of spacelike momenta. We 
shall concentrate on class IV contributions exclusively as we are interested in computing 
the impact of the new isospin-violating sources on AA. The knowledge of the q^ = 
couplings now no longer suffices. One is forced to model the momentum dependence of the 
coupling constants — including that of the isospin-violating components. Here we consider 
two different estimates for the dependence of the CSB potentials. First, we simply adopt 
the momentum dependence which emerges from fits to the isospin-conserving two-nucleon 
data. Thus, the ratio of the isospin-violating to the isospin-conserving coupling, e.g., fi /g^, 
remains unchanged. Note that in the Bonn model fi/gl is also a constant. We implement 
this choice by modifying the meson-nucleon "point" couplings indicated in Eq. (0) as per 
the Bonn B potential parameters, see Table |T|. That is. 
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9nnp 



(28a) 
(28b) 



This is an additional model assumption. Here we use q to denote the three-momentum 
transfer; we consider the form factors in the Breit frame, where qo = and = — q^. 
Second, we compute the 0(q^) isospin-breaking in the couphngs in the nonrelativistic quark 
model, in order to gauge the uncertainty in the momentum dependence of the CSB potentials. 
Let us examine the isospin-breaking in the Sachs- Walecka form factors separated into 



contributions from the isoscalar or isovector quark charges. These quantities are related to 
the u and p couplings by virtue of our vector dominance assumption. As previously, we will 
discuss merely the isospin breaking in the NNu vertex in detail. Now 



G 
G 



E,p 
E,n 



'^9d{d)n + 9u{u)n , 



and 



G' 



M,p 



2M„ 



M,n 



2M„. 



18 \m„ / 

^ (—{d)n - —{u)„) . 
18 \md rriu / 



(29a) 
(29b) 



(30a) 
(30b) 



These expressions are generalizations of Eqs. (P) and (0). We have used the notation of 
Eq. (ID in denoting the isoscalar and isovector quark charges and have introduced {u)p, for 
example, to represent the Fourier transform of the proton wave function with respect to 
the up quark coordinate. We compute the latter in the harmonic oscillator quark model 



for simplicity. In the harmonic oscillator quark model ||21| , ^ the nucleon possesses a mass 
Mjsi = 2mi +m2, so that for the proton mi = m„ and m2 = m^. For convenience one defines 

2_ foizzT' „,i ™ — 317111712/ {2mi + 1712), mp 



Sknip and R^^ 



y/3kmx, where mx 
spring constant. One finds that 



mi, and k is the 



{u)p = (exp(zq- r„))p 

{d)p = (exp(zq- rd))p 
{u)n = (exp(«q- r„))„ 
{d)n = (exp(zq • rd))„ 



1-3! 



3q" 



Rpp + 3 



m 



1 - 



q 



md\ 2 
.Mp) ""'^ 



2 / N 2 



2 / 



+ 0(q^ 



(31a) 

(31b) 
(31c) 
(31d) 



We write the Fourier transforms in Eq. ( pi]) through O(q^) only. This suffices to make 
contact with the hadronic form factors. Moreover, one cannot expect the nonrelativistic 
quark model to be reliable at still larger momentum transfers. We must now relate the 
above electric and magnetic form factors to the /'s and g^s present in the definition of the 
vertex, Eq. (|l]). Following the usual relation between the electromagnetic form factors Ge, 
Gm and Fi, F2, vector dominance dictates that 
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(32a) 
(32b) 



Fits to the electromagnetic form factor data indicate that -Fi(g^) and -^2(0'^) fall with different 



rates in g^; vector dominance implies that this should be true of g'^^ and f'j^ 



as well. Note, 



this is at odds with the Bonn model, as it assumes that the ratio / gti'' is constant. 
We proceed as follows. We compute ffj and gfj to (9(q'^, Am^), using Eqs. (p^-|3~ 



\ ^^) . Then 

we estimate the "effective" A^, as defined in Eq. (p^), required to reproduce the isospin 
breaking computed to 0(q^, Am^) and use that A^^ in our subsequent computation of the 
spin-singlet-triplet mixing angles. Thus, 



9^ + 9ir, 







Am 

+ 

m 


" 5q2 




6 


24M2 


5 Am 


1 - 


q2i?2 


1 


6 m 




3 


4M2 



q2i?2 



72 

+ 0(q^ Am^ 



+ 0(q^ Am' 



(33a) 
(33b) 



Several remarks are in order. First, note that we have defined = V^ikm, where m is the 
average mass of the up and down quarks. From Eq. ( |33b| ) we observe = to C(q*^, Am^); 
this is consistent with the Bonn model, which assumes = for all q^. We have performed 
the same calculations for the A^A^p vertex as well. In this case, one finds results at odds 
with the Bonn model, as fi/gl is not constant to 0{c^ , ^m^). Note that at nonzero q' 
CSB potentials beyond those enumerated in Eqs. (|19a|) - (|19d| ) may exist. For example, at 
(9(q2,Am) a new CSB contribution arises from the combination fig^. Yet, like the rho 
contribution to the CSB potential given in Eq. (|19b| ), it is not numerically important, due 
to the small value of fi'jvAfp Bonn model — recall that fi'^^^/S'^^^ ~ 27 |2^. Let 

us proceed to examine the impact of Eq. (|33b| ) on the omega contribution to the class IV 
CSB potential. We fix the scale R by requiring that the isospin conserving vertex, g^, fall 



in q at the rate given by the Bonn model, so that R = \/l2/A^ ^ .37 fm. We choose R 
in this manner as our primary interest is in determining the fall-off of the isospin-breaking 
potential relative to the isospin-conserving one. Noting Eq. ( |19a|) , we consider 



fi9^ 



5 Am 



6 m 
5 Am 



1-6 



q^ 



4M2 / 



+ C>(q^ Am' 



6 m 



(34a) 
(34b) 



By replacing the A^^ of Eq. ( pHf ) with the A^^ given above, such that 

A'. = A' ^ ^ ^ 



^6 + A2/4M2 



(35) 



we obtain an expression for the CSB potential, Eq. (|19a|) , which is of the form given by our 
original prescription, Eq. (^), yet is equivalent to the isospin-breaking calculated in the 
harmonic oscillator quark model at (9(q2,Am). Numerically, the Bonn model A^^ = 1850 
MeV is changed to A^^ = 1401 MeV. At this order the coefficient of g1 is not negative, so 
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that we cannot carry out the above exercise for the rho as welL The rho's numerical impact 
on the CSB potential is small, so that this gap does not impact our uncertainty estimate in 
any significant way. We will proceed to compute the spin-singlet -triplet mixing angles for 
the potential given by Eqs. ( |19a|) and (pSD for both the A^^ of the Bonn potential and the 



A,; of Eq. 

In Fig. we present estimates of the CSB potentials given in Eqs. (|19a|) , (|19b| ), and (|21a]) 



using Eq. with the A^^^p of the Bonn model. The qualitative conclusions we draw here are 
not sensitive to the choice of A^ p, so that we simply present the potentials computed in the 
Bonn model. The solid line is the CSB potential which results from p-uj mixing, Eq. ( ^la| ), if 
the on-shell value of the p-uj mixing amplitude is employed for the entire range of momenta. 
This is the potential traditionally used in studies of CSB observables. A potential of this 
strength is required to describe the analyzing power difference A A measured in elastic n — p 
scattering |jl6|]. In contrast, the dashed line results if the momentum-dependent p-u mixing 
amplitude of Ref. is employed in Eq. ( |21a| ) — this is too small to fit the data [0, yet a 
model in which the vector mesons couple to conserved currents must yield a vanishing mixing 
amplitude at = . We have not extended our model to describe p-u mixing; the vector 
dominance assumption we use implies, however, that the = mixing must be zero in 
this framework. We take the momentum- dependence of the mixing amplitude computed by 
Piekarewicz and Williams as archetypal. This latter CSB potential in itself would upset 
the previous agreement with experiment. However, the new sources of isospin violation 
computed here are sufficient to restore the agreement. In particular, the contribution from 
omega- meson exchange, given by the dash-dotted line, is large and comparable in magnitude 
to the one arising from on-shell p-u mixing. We have also computed the contribution from 
the rho-meson, given by the dotted line, though it is negligible due to the small NNp vector 
coupling. 

In Fig. ^ we display the above CSB potentials in configuration space. The potentials 
have been normalized so that the areas under the curves equal Vlq"^ = 0). Qualitatively, the 
trends observed in Fig. ^ remain: we obtain large contributions from on-shell p-u mixing 
and omega-meson exchange and small corrections to the latter from off-shell p-u mixing and 
rho-meson exchange. These results are suggestive, yet we can obtain a precise estimate of the 
impact of the enumerated isospin-violating sources on AA by calculating the spin-singlet- 
triplet mixing angles, 7j. These are the dynamical quantities driving AA |TO| , p!T| . Recall 
that the elastic scattering amplitude of two spin-1/2 particles is specified by six invariant 
amplitudes a, b, c, d, e, and / [0], so that 

1 r 



M 

2 



a + b) + {a - 6)((Ti ■ n)(cr2 • n) 

+ {c + d) ((Ti • m) ((72 ■ m) + (c - rf) (cri ■ \) (^2 ■ 1) (36) 
+ e(cTi + a2) ■ n + /(ai - 0-2) ■ n , 

1-^^. "'^wi^r "^w^i- 

and kj and kf are the initial and final cm. momenta of particle 1. The -n — p analyzing 
power difference is nonzero only if accompanied by spin-singlet-triplet mixing, specifically 



where 
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AA(^) = MO) - AM = 2Re{b*f)/ao , (38) 

where (Tq is the unpolarized differential cross section. The spin-singlet-triplet mixing is 
controlled by /. Neglecting electromagnetic effects, / is connected to the mixing angles 7j 
via 
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■ oo 



/(^' ^) = ^ E (2^ + 1) M^lj) eM^Sj + i5jj)dUe) , (39) 



J=l 



where the df^lO) are Wigner functions and the 6j and 6jj are the singlet and uncoupled 
triplet bar phase shifts, respectively. In a distorted-wave Born approximation the mixing 
angles themselves are given by flOf 



y roo roc 

J J = -mk^J{J + l) drr^Rj{r)Viv{r)Rjj{r) = drlj{r) , (40) 
where we have introduced the class IV CSB potential 

Note that Rj{r) and Rjj{r) are the spin-singlet and triplet radial wave functions, respec- 
tively, for NN scattering in the L = J channel. Distortion effects are incorporated through 
these radial wave functions; we assume them adequately described by solutions to the Reid 
soft-core potential |^ . In Table |^ the first four nonvanishing mixing angles, J = 1 — 4, are 



presented at a laboratory energy of 183 MeV. In addition, the integrand from which 71 is 
obtained, that is, /i(r) in Eq. (|40|) , is plotted in Fig. |^. This represents the class IV potential 
suitably weighted by realistic A^A^ wave functions. Three calculations are presented for com- 
parison. The solid line is obtained using Eq. (|21a|) and the on-shell value of the p-u mixing 



amplitude; the area under this curve is the mixing angle required to reproduce the AA data. 
In the dashed line we have combined the off-shell p-u mixing contribution described above 
with the isospin- violating vertex contributions arising from omega- and rho-meson exchange. 
Albeit form factor uncertainties in the isospin-violating vertices, this is our best estimate of 
the mixing angle contribution. The vertices in this figure were evaluated using Eq. (^) and 
the Bonn cutoff parameters A^^ p tabulated in Table |IT1[ We have also combined the on-shell 
p-uj mixing contribution with the above vertex contributions, even if our model is not con- 
sistent with a nonzero mixing amplitude at = — this is shown by the dash-dotted line. 
The integrand in this case is considerably larger than the other two estimates. The J = 1 
mixing angles for these integrands are displayed in parentheses next to the curve labels. 
The agreement between the first two calculations is very good. Indeed, the contribution 
from omega-meson exchange, together with small corrections from off-shell p-uj mixing and 
rho-meson exchange, results in a 3% reduction in the value of 71, relative to the on-shell 
value. This kind of agreement — at the few percent level — is maintained throughout all 
the examined partial waves, note Table |IV|. These computations have also been performed 
with the form factor A^, Eq. (^), estimated in the harmonic oscillator quark model. The 
mixing angles obtained in this fashion vary by about 10% in the important partial waves 
from those computed with the Bonn form factors; note that 71 = .036°, rather than .033°. 
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For a detailed comparison, see Table |^ — the mixing angles which use the harmonic oscil- 
lator quark model results to estimate the dependence of the CSB potentials are shown 
in parentheses. The mixing angles computed with isospin-breaking meson-nucleon vertices 
and off-shell p-uo mixing in the two approaches bracket the old on-shell p-uo mixing results 
for J = 1 — 3, so that these new estimates are also quite close to the "old" on-shell results. 
Ku) is some 3/4 of A^, yet the above calculations show that the at 183 MeV is essentially 
dominated by the = physics. Note that if one were to assume a momentum-independent 
p-uj mixing amplitude |jl3[ and to include the contributions from omega and rho exchange 
an increase of almost a factor of two relative to the above mixing angle estimates would 
result. 



V. CONCLUSIONS 

We have studied the charge-symmetry breaking in the NN potential arising from isospin- 
violating meson-nucleon coupling constants. The isospin-violating couplings are obtained by 
computing matrix elements of quark currents of the appropriate Lorentz and flavor structure 
between nucleon states. We have used a nonrelativistic quark model to evaluate these matrix 
elements, yet our estimates at = depend merely on the spin and flavor structure of the 
nucleon wave function, rather than on the details of the quark momentum distribution. Thus, 
in the vector meson sector, for example, our model estimates at = depend on our vector 
dominance assumption, but little else. We have also studied isospin breaking in the NNn 
and NNa vertices. No isospin- violations exist in the a vertex at = 0. We have found that 
the breaking in the NNtc vertex depends on whether the ttN coupling is of pseudoscalar or 
pseudovector character — no isospin breaking results if pseudovector coupling is assumed. 
However, a pseudoscalar vrN coupling is commonly used in studies of CSB, and the breaking 
we find in the vertex is substantially larger than the breaking computed in hadronic models 
of neutral pion exchange. Thus, any CSB observable receiving an important contribution 
from TT-T] mixing will also be affected by the exchange of neutral pions. 

We have found that omega-meson exchange is an important component of the class IV 
charge-symmetry-breaking NN potential needed to describe the analyzing power difference 
measured in elastic n — p scattering at low energies. The potential which emerges from 
the isospin-violating NNu vertex is identical in structure to that from p-u mixing |]15| . 



Moreover, our estimates indicate that these two contributions — with the mixing amplitude 
fixed at its on-shell value — are comparable in magnitude and identical in sign at g2 = 0. 
Models in which the vector mesons couple to conserved currents, of which ours is an example, 
have no p-u mixing at g^ = [Q. We have found that isospin- violation in the NNu 
vertex can generate a CSB potential of sufficient magnitude to fill the phenomenological 
role required by the lUCF measurement of A A at 183 MeV. 

The isospin-violating couplings we have computed at g^ = do not suffice to make a 
quantitative prediction of the CSB potential needed for the lUCF experiment. One must 
compute the g^ dependence of the NN meson vertex as well — including that of the isospin- 
violating pieces. We have considered two simple estimates. The first is simply a prescription: 
we modify the "point" couplings by introducing hadronic form factors according to the Bonn 
B potential. This assumes that the relative strength of the isospin-breaking potential found 
at g^ = persists at finite g^ as well. In the second we compute the isospin breaking to 
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(9(q^, Am^) using the spatial wave functions of the harmonic oscillator quark model and 
find the hadronic form factor for omega exchange needed to reproduce the isospin breaking 
computed to the above order. The use of the spatial component of the nucleon wave function 
is required here, so that this estimate is rather more model dependent than our — Q results. 
We find that the use of the latter estimate yields slightly larger CSB potentials. 

Armed with estimates of the momentum dependence of the NN meson vertex, we have 
computed the spin-singlet-triplet mixing angles 7^: these are the fundamental dynamical 
quantities driving ^A. Our 7^ computation is realistic as we have used the Reid soft-core 
potential to generate the distortions in the NN weive function. We have computed the 
spin-singlet-triplet angles using three different sources of isospin violation: (1) p-cu mixing 
with the amplitude fixed at its on-shell value, (2) off-shell p-cu mixing plus omega- and rho- 
meson exchange, and (3) on-shell p-u mixing plus omega- and rho-meson exchange. The 
first case, used in the original estimates of AA, represents a "baseline" value, as it fits the 
data. A CSB potential of this magnitude accounts for a substantial fraction of the measured 
value of AA at 183 MeV. The second case, which should be regarded as our best estimate, 
yields values for 7j that are within 10% of those obtained with on-shell p-oo mixing, for the 
important partial waves. In contrast, case (3) results in a factor-of-two enhancement relative 
to the original calculation using on-shell p-uo mixing. Two important results thus emerge 
from the present work. First, we have found a new source of isospin violation, namely in 
the NNlo vertex, which can fill the role demanded by the data. Second, we have shown that 
insisting upon a p-uj mixing amplitude held constant at its on-shell value, after including 
the contribution from omega-meson exchange, results in a class IV potential too large to be 
consistent with the lUCF AA measurement. 
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FIGURES 



FIG. 1. Charge-symmetry-breaking component of the NN potential as a function of arising 
from on-shcU p-to mixing (sohd hne), off-shell p-uj mixing (dashed line), omega- meson exchange 
(dash-dotted line), and rho-meson exchange (dotted line). 

FIG. 2. Charge-symmetry-breaking component of the A^A^ potential as a function of r arising 
from on-shell p-u! mixing (solid line), off-shell p-cu mixing (dashed line), omega- meson exchange 
(dash-dotted line), and rho-meson exchange (dotted line). 

FIG. 3. The integrand of the spin-singlet-triplet mixing angle 71 at 183 MeV for three different 
estimates of CSB: on-shell p-u mixing (solid line), off-shell p-oj mixing plus omega- and rho-meson 

exchange (dashed line) , and on-shell p-uj mixing plus omega- and rho-meson exchange (dash-dotted 
line). The value of 71 for each estimate, which is simply the area under the appropriate curve, 
appears in parentheses next to its label. 
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TABLES 



TABLE I. Quark model charges. The superscripts denote scalar, pseudoscalar, vector-isoscalar, 
and vector-isovector quark charges, respectively. 
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9' 




5+ 


9~ 
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+1/3 


+1/5 




+1/3 


+1 
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+1/3 
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+1/3 


-1 


TABLE IL Proton, neutron, isoscalar, and isovector 


meson- 


■nucleon coupling constants. 
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TABLE m. 


Meson masses, coupling 


: constants, tensor-to-vector ratio. 


and cutoff parameters 


of the Bonn B potential. 










Meson 


Mass (MeV) 






f/9 


A (MeV) 


TT 


138 


14.21 






1700 


V 


549 


2.25 






1500 


p 


769 


0.42 




6.1 


1850 




783 


11.13 




0.0 


1850 



TABLE IV. Spin singlet-triplet mixing angles 7j (in degrees) at a laboratory energy of 
Tiab = 183 MeV. The values in parenthesis use a form factor computed from the quark model 
(see text for details). 



J 


{p\H\u)\on 


{p\H\uj)\oS + {u;- p) 


{p\H\u)\on + {UJ- p) 


1 


3.41 X 10-2 


3.32 X 10-2 (3.66 x 10-^) 


6.70 X 10-2 (7.22 X 10-2) 


2 


4.51 X 10-2 


4.40 X 10-2 (4.88 X 10^2) 


8.77 X 10-2 (9.52 x 10^2) 


3 


3.77 X 10-3 


3.71 X 10-3 (4.71 X 10-3) 


6.00 X 10-3 (7.55 X 10-3) 


4 


8.04 X 10-^ 


8.58 X 10-^ (1.12 X 10-3) 


1.18 X 10-3 (1.55 X 10-3) 
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